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ABSTRACT 

Spatial maps of the high-pass frequency filtered time-averaged root-mean-squared (RMS) Doppler velocities tend to 
show substantial decrements within regions of strong field and curiously, randomly distributed patches of enhancement 
in the vicinity. We propose that these haloes or enhancements are a consequence of magnetic-field-induced mode mixing 
(scattering), resulting in the preferential powering of waves that possess strong surface velocity signatures (i.e. scattering 
from low to high wavenumbers) . Evidently, this process can occur in the reverse, and therefore in order to determine if 
the haloes are indeed caused by mode mixing, we must answer the question: how are acoustic waves scattered by magnetic 
fields? Through simulations of the interactions between waves and sunspots and models of plage, we demonstrate that 
the high to low modal order scattering channels are favoured. With increasing frequency and consequently, decreasing 
wavelength, a growing number of modes are scattered by the sunspot, thereby rendering the enhancements most visible 
around the high-frequency parts of the spectrum. The haloes obtained from the simulations are on the same order of 
magnitude but weaker than those observed. We also present observational evidence to support this theory: observations 
of active region AR9787 are firstly frequency filtered to isolate the 5-6 mHz signals and secondly, decomposed into three 
wavenumber bandpasses, I — [0, 400], [400, 800], [800, 2222]. With increasing wavenumber, the extent of the halo effect is 
seen increase dramatically, in line with theoretical expectation. 
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1. Introduction 

Observations of acoustic 
netic field are numerous 
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ower in re g ions o f strong mag- 
Lites et all (Il982h: iBrown et al.1 
Bogdan et all dl993h 
Hindman fc Brownf (| 199S 



Donea e tall (l2QQ0h: I Jain fc Haberl (|2QQ2h : IMoretti et all 
(|2007t ): iNagashima et al.l (|2007l ). In general there is agree- 
■ ment that trapped acoustic wave power (with frequencies 
k/H below the acoustic cut-off) measured using Doppler veloc- 
ities is strongly suppressed in the strongest field regions, 
d ' and that in the weak field surrounding area, there is an 
enhancement of wave power at frequencies comparable to 
and above the acoustic cut-off. Corresponding measure- 
ments in intensity however do not quite provide such a 
consistent picture; numerous differences are seen in the 
structure and spatial locations of the haloes. The variety 
of instruments and optical spectral lines used in taking 
these observations may b e the source of the discrepancies 
(e.g. IMoretti et all l2QQ7h . Furthermore, because intensity 
is a complex metric, and interpreting changes thereof can 
be non-trivial, we stick to the Doppler measurements. 

There have been several theorie s attempting to ex - 
plain the emergen ce of the haloes. [Brown et al.l ((1992); 
iDonea et al.l (|2000l ) subscribe to the hypothesis of excess 
wave source activity in the vicinity of active regions. Such 
an assumption is entirely plausible given the remarkable as- 
sociated changes in th e turbulent laye r s of t he solar photo- 
sphere. More recently, Ijacoutot et al.l ([2008), based on nu- 
merical simulations of magnetoconvection found that with 



increasing strength of network magnetic fields, the oscil- 
lation spectrum shift ed t o higher frequenc ies. Alternately, 
IKuridze et all (|2Q08h and [H anasogd (|2008l ) have proposed 
that a linear (i.e. small wave amplitudes) magnetohydro- 
dynamic (MHD) wave interact ion mechanism m ay be re- 
sponsible for the haloes. While iHanasogej ([2008), based on 
numerical simulations of MHD wave propagation, merely 
speculated about t he existence of such a mechanism, 
IKuridze et al.l (|2008f ) suggested that overlying canopy mag- 
netic fields may locally raise the acoustic cut-off frequency, 
thereby trapping the hi gh frequency wave s . On a related 
note, the observations of iBraun fc Lindsevl (|2QQQh indicate 
that high frequency waves undergo enhanced reflections in 
active regions. 

Continuing on the path set forth in lHanasog 1 (120081 ). 
we posit here that wave scattering induces a transfer 
of energies across modes at constant frequency, a conse- 
quence of which is the appearance of the velocity ampli- 
tude enha ncements in the v icinity of active regions. Very 
recently, iGordovskvv et all (|2009f h in a study of p-mode 
oscillation amplitude anomalies around active regions (e.g. 
iHindman et al.| [l997) reached a similar conclusion, i.e., that 
the scattering occurs preferentially from high to low radial 
orders. In this contribution, we dynamically solve the lin- 
earized three-dimensional MHD equations in order to ex- 
amine the mode interactions with strong-field flux tubes 
and plage. The principal assumption in these calculations 
is that of linearity, or low wave amplitudes in comparison 
to the local phase and group velocities. Analyses of the 
simulated data display signs of biased scattering and the 
consequent emergence of velocity enhancements. We also 
present observational evidence to support this theory. 



2 



Hanasoge: Haloes by Wave Scattering 







1000 f 



100 r- 




[mHz] 



4 



1 



500 



1000 

I 



1500 



2000 



Fig. 1. Arbitrarily scaled mode mass (M n ) as a function 
of frequency for different p n ridges compute d for a poly- 
trope (equation (3.13) of Bogda n et al.| [l996). The curves 
from bottom to top correspond to the n = 1, .., 5 p n ridges. 
The ratios of the masses are useful indicators of the de- 
gree of change in the surface velocity signatures that can 
be effected by the scatter from one mode to another. 



2. Theory 

A wave propagating in the quiet Sun has an associated dis- 
placement (for example) eigenfunction. When this wave en- 
counters an anomaly, the wave displacements may no longer 
be well represented by just the quiet Sun eigenfunction, but 
rather, require an additional basis of eigenfunctions. These 
additional eigenfunctions can only be introduced into the 
system by the creation of a corresponding set of waves. 
The energy required to produce these additional or scat- 
tered waves is extracted from the incident wave, resulting 
in an energy loss of the original wave. Due to the low ampli- 
tudes of the seismic waves of relevance to helioseismology, 
we may treat scattering in the Sun as a linear phenomenon. 
Furthermore, the background medium and associated scat- 
terers are considered to be time stationary. Consequently, 
wave interactions with scatterers can only result in the re- 
distribution of energies across wavenumbers at constant fre- 
quency. 

The velocities observed at the solar surface are a sum of 
the velocities of a large number of modes, each of which has 
a characteristic surface amplitude. The ratio of the time- 
average energy of a mode to the surface mean squar e ve- 
locity is termed mode mass (e.g. Bo gdan et al.l [l996h : low 
mode mass waves possess strong surface signatures and con- 
versely. Mode mass is a monotonically decreasing function 
of wavenumber and frequency, as seen in Figure [TJ Consider, 
for example, a p^ mode of unit energy being fully trans- 
formed into pi at a specific frequency. From Figure 1, we 
see that the p4 to p\ mode mass ratio is approximately 
5, indicating the magnitude by which the surface velocity 
square would increase were such a transformation to occur. 
A small amount of scatter can thus result in significant 
changes in the observed surface Doppler velocities. 

Observations of the quiet Sun photosphere show a 
modal velocity spectrum that peaks around a frequency 
of v ^ 3.6 mHz. This is the consequence of an impasse 
reached between damping, emission, and mode mass. The 



Fig. 2. Log power spectrum from 1440 minutes of full disk 
observations of the active region AR9 787 by the Michelsq n 
Doppler Imager instrument (MDI; IScherrer et al.l Il995l ). 
The frequency is denoted by v and the spherical harmonic 
degree by I = 27rR@/\, where A is the wavelength. The 
vertical line marks 2ttRq/I c ^ where l c = 20 Mm, the char- 
acteristic length scale of a hypothetical scatterer. We may 
broadly classify the scattering process into two regimes, 
weak and strong, on the left and right sides of the line 
respectively. Roughly speaking, modes on the left of the 
line do not "see" the scatterer very well because of their 
relatively-large wavelengths. At high frequencies, the wave- 
length decreases and many more modes couple with the 
scatterer, potentially significantly altering the distribution 
of energies with wavenumber. Note that at 3 mHz, four 
ridges lie in the "strong" scattering regime whereas at 5 
mHz, 11 radial orders fall in this regime. 



presence of strong scatterers can disturb this delicate equi- 
librium and reconfigure the spectrum. This, we propose, is 
precisely what is happening in the vicinity of active regions. 
High mode mass waves are depleted of their energies, which 
are then channeled into the low mode mass counterparts. 

Consider a hypothetical situation where waves are scat- 
tered by a flux tube of size 15 Mm or so. Waves with wave- 
lengths much larger than this flux tube do not couple with 
it very strongly and may be classified as being weakly scat- 
tered. Conversely, smaller wavelengths (also higher frequen- 
cies) are strongly scattered. It is therefore reasonable to 
conclude that with increasing frequency, more modes fall 
into the strong scattering regime. This leads to the modal 
energy being redistributed to greater extents at higher fre- 
quencies (Figure EJ). Above the acoustic cut-off frequency, 
however, the drop in acoustic power makes observing haloes 
more difficult, requiring the use of other spectral lines or 
high cadence/resolution observations. By choosing Sodium 
and Potassium li nes that form 500 an d 250 km above the 
photosphere, e.g. iMoretti et al.l (|2007l ) have observed the 
appearance of haloes at frequencies as h igh as 7 mHz. 

Hankel analyses of observati ons (e.g. Braun et al.lll992l : 
Bog dan et al.|[l993l ; lBraun|[l995l ) contain all the features we 
might expect from this mode mixing theory: (1) the phase 
shifts increase strongly with frequency, (2) the absorption 
coefficients are seen to initially increase with frequency and 
saturate at high frequencies for the ps and higher ridges 
while for the low mode mass ridges (^1,^2), the absorption 
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falls to zero at the high frequencies. Naively, we might con- 
clude that the reason for this drop in the absorption is an 
excess of low mode mass waves introduced by scattering 
from high mode mass ones. However, it turns out that the 
damping rates of the pi,P2 modes at high frequency are 
so large that the analyses are dominated by the local ex- 
citation and decay of the modes, as opposed arising from 
correlations across the sunspot (private co mmunication, A. 
C. B irch 2009; also, Figures 13 and 14 of lBraun fc Birchl 
2008). We are therefore unable to conclusively state a case 
for this theory from just the Hankel observations. 

Thus far, we have only discussed the case when the high 
mode mass waves scatter into the lower ones. A double- 
edged sword now awaits us: scattering is a two-way street. 
Equally likely to happen is scattering from low mode mass 
waves into the high, resulting in a surface velocity deficit. 
The question then is: does one form of scattering (high to 
low) preferentially occur over the other? In order to answer 
this, we compute the linear interaction of waves with (a) a 
sunspot sized flux tube, and (b) a bipolar plage-like cluster 
of moderately sized flux tubes. By ridge filtering the simu- 
lated data, we separate the contributions of the different p 
ridges to the halo. We do not present corresponding anal- 
yses of the observations because of the low signal-to-noise 
properties. 

The three-dimensional ideal MHD equations are solved 
via the procedure described in iHanasogej (|2008l ). For the 
sunspot calculation, a flux tube with peak photospheric 
field strength of 3000G and radius of 10 Mm (full width 
at half maximum) at the photosphere is placed at the hori- 
zontal centre of a box of size 200 x 200 x 35 Mm 3 (horizon- 
tal x horizontal x depth) resolved using 256 x 256 x 300 
grid points. The d etails of the setup are described in 
iMoradi et al.l (|2009l ) . The system is deterministically forced 
with dipolar vertical velocity sources. The sources are sup- 
pressed in the sunspot region so as to avoid exciting waves 
in regions of strong magnetic field. In order to determine the 
extent of the halo effect, we compare the MHD calculation 
with a quiet simulation, the latter corresponding to a com- 
putation of the wave propagation in a domain with sources 
everywhere and no sunspot. The plage calculation consists 
of five moderate to weak bipolar flux tubes placed at ran- 
dom locations around the horizontal centre of the simula- 
tion box. The plage model used in these calculations fails at 
mimicking the spaghetti-like complexity of solar plage and 
we may therefore b e significantly underestim ating the de- 
gree of scatter (e.g. lHanasoge Callyl l2QQ9h . Because the 
plasma- /3 at the source excitation depth is much greater 
than 1 in the plage case, the sources are chosen to be as in 
the quiet simulation, i.e., uniform everywhere. In Figure OS 
we display the magnetic field distributions of the sunspot 
and plage models at a height of 200 km above the photo- 
sphere. Vertical velocity components are extracted at an al- 
titude of approximately 200 km above the photosphere (the 
upper turning point of the waves in the quiet simulation is 
the photosphere). This is done so as to mimic the quiet Sun 
height of formation of the Nickel line that MDI observes. 
However, the signals are sampled at this geometric height 
in all regions including in areas of strong magnetic field. 
We are careful therefore not to interpret the signals in the 
regions where the field may be significant. It is also impor- 
tant to note that we use a stratification quite different from 
model S, and consequently the calculations do not contain 
the same resonant modes as the Sun. 
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Fig. 3. The magnetic field distributions of the sunspot (left 
panel) and the bipolar 'plage' models (right). The dots 
mark the location of the half maximum of the field strength 
in the sunspot model. 



Three caveats must be kept in mind when interpret- 
ing the results of these simulations: (1) for the background 
stratification, we use a combination of sub-surface poly- 
tropic and atmospheric isothermal layers, meaning that the 
resonant frequencies differ from the solar values, (2) we do 
not incorporate realistic mode damping, which means the 
waves are correlated over larger distances, and (3) an Alfven 
speed limiter is applied in order to prevent the plasma- in 
the atmosphere from falling too low (e.g. lHanasoge! [20081 ) 
- this avoids negative pressures and lifts severe timestep 
restrictions. 

We pause the discussion for a moment now to define 
the terms absorption and scattering. The former describes 
a process by which acoustic waves undergo transformation 
into slow or Alfven MHD waves (e.g. lCallv fc Bogdan|[l997[ ) 
and leak out into the atmosphere or disappear into the in- 
terior. They are not a part of the resonant ridges of the 
observed power spectrum. The latter denotes mode mix- 
ing between the various acoustic modes. In this context, 
the Alfven speed limiter described above is a cause for 
some worry since it can alter the balance between absorp- 
tion and scattering. Through preliminary work with Hankel 
transforms (private communication, S. Couvidat, 2008), we 
have found that the sunspot in our simulations acts more 
as a scatterer and less of an absorber. This is to be con- 
trasted with the work of ICameron et all ([2008), who find 
that the sunspots in their simulations absorb very effec- 
tively. However the calculations with plage are more be- 
lievable because the plasma- (3 in the tubes become very low 
only at regions close to the upper boundary. Consequently, 
the Alfven speed limiter acts only near the boundary and 
possibly does not affect the solution within. 



3. Analyses and Conclusions 

The power spectrum of the simulated data is shown in 
Figure [4^. The high-frequency filtered power maps of 
the two simulatio ns (panels b,c) and t he observations o f 
AR9787 (panel d; IScherrer et all 119951 : iGizon et al J l2QQ9h 
show visible haloes. In order to reduce the noise in the 
analyses of the simulated data, we also compute a quiet ver- 
sion, i.e. using the same realization of the source function 
but with no sunspot and sources everywhere. The filtered 
power maps of the quiet cube are then subtracted from the 
MHD counterparts (see appendix [Al for a description of the 
filters). The frequency filtered maps of the simulations in 
Figured] (b,c) reproduce some of the observational trends - 
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Fig. 4. Log modal simulation spectrum (panel a) and power 
maps of the 5-6 mHz frequency-filtered simulated sunspot 
(panel b), bipolar 'plage' (panel c), and full-disk observa- 
tional data (panel d); the scales are saturated in all the 
plots. We do not force the system much below the / ridge, 
a manifestation of which is the artificial drop in power in 
panel a. In panels b and c, noise subtraction is applied in 
order to gen erate a clear view of the ac oustic power (e.g . 
IWerne et all 12004 lHanasoge et~aD I2QQ7I : E anasog e 2008). 
The quiet calculation contains no magnetic field and has 
a uniform source distribution over the entire horizontal ex- 
tent of the domain (i.e., no source suppression in the strong 
field regions). 



amplitude depressions in the strong field regions and haloes 
at high frequencies. 

Shown in Figure [5] are the power maps and flux-tube- 
centred azimuthal averages of the high-pass- andpi,p2, and 
Ps ridge-filtered data. The effect described in the previous 
section is seen quite clearly - namely that the p\ , P2 ridge- 
filtered maps show haloes while a power decrease persists in 
the ps map. Strong scattering has caused high mode mass 
power to be channeled into the lower mode mass waves. 

The most promising demonstration of the validity of 
this theory is displayed in Figure [6l We take the high- 
frequency filtered observations of AR9787 and decom- 
pose it further into three different spatial bandpasses: 
I - [0,400], [400,800], [800,2222], where I is the spherical 
harmonic degree (see appendix O for a description of the fil- 
ters). Theoretically, we expect to see increasing halo power 
with spatial bandpass, since at a given frequency, the low- 
est mode mass waves are at the high wavenumbers. As the 
quiet Sun control region, we choose a ~ 147 x 61.3 Mm 2 
sized area, essentially the upper left part of Figure Sfd). 
The power in each bandpass is normalized by the average 
power value of the correspondingly filtered quiet Sun region. 
We then compute the areal power distribution: the number 
of pixels that contain a specific value of relative power are 
counted and a histogram is plotted. Theoretically, we ex- 
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Fig. 5. High-pass (5-6 mHz), pi,P2,P3 ridge-filtered power 
maps of the linear MHD simulations of waves interacting 
with a sunspot sized flux tube. A halo is seen around the 
pi ridge, the p2 map is flooded with excess power outside 
of the strong field regions, and the p% ridge shows signifi- 
cant power depletion throughout. The panels on the right 
show the azimuthal averages around the centre of the flux 
tube. The y range of the azimuthally averaged p\ power 
plot has been restricted. Identically filtered power maps of 
a quiet simulation (i.e. with no sunspot and sources every- 
wher e) have been subtracted in order to r educe the noise 
(e.g. lHanasoge et al.ll2007t I Hanaso gel 120081 ). 



pect that with growing wavenumber, the enhancement area 
should increase. And indeed, this is exactly what is seen in 
Figure Eta)! The number of pixels showing an enhancement 
greater than 1-<j, where a is the variance of the power map, 
in the three spatial bandpasses is (0.25,0.31,0.38) x 10 6 
respectively. We also show the quiet Sun and active re- 
gion power distributions of two bandpasses ([400, 800] and 
[800,2222]) in panel b of Figure [6] in order to demonstrate 
that with increasing wavenumber, the distribution shifts in 
the direction of increasing power. 

The theory presented here is capable of explaining the 
observations of haloes around active regions. The combina- 
tion of a sunspot and surrounding plage can significantly 
alter the acoustic wave power in the vicinity, thereby shift- 
ing the local oscillation power spectral density. The magni- 
tudes of the enhancements seen in the simulations are seen 
to be not dissimilar from the observationally derived val- 
ues - this implies that the mechanism proposed here, with 
suitable improvements, may be able to reproduce the mag- 
nitudes of the observed haloes. 
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Fig. 6. High-pass (5-6 mHz) and wavenumber filtered ob- 
servations of AR 9787: shown in panel a is the number of 
pixels in the active region power map over which a given 
relative power level persists; in b, we plot the power dis- 
tributions over the quiet Sun (QS) and active region (AR). 
The power in each bandpass is normalized with respect to 
the mean over a correspondingly filtered quiet region; here 
the quiet Sun control is chosen to be a 147 x 61 Mm 2 sized 
region on the upper left part of Figure ffl(d). With increas- 
ing wavenumber, larger numbers of pixels display enhance- 
ments. The power distribution in (b) is also seen to grow 
progressively asymmetric; the red curves show the active 
region power distribution in the two bandpasses. 



The analyses of the observations show that the high- 
wavenumber, high-frequency parts of the spectrum are the 
most affected, containing the strongest velocity enhance- 
ments. The trend that we sought, and recovered, was that 
of progressively increasing halo amplitudes with wavenum- 
ber. With more sophisticated analyses, we may be able to 
extract information about mode mixing and perhaps refine 
the interpretations of hankel absorption coefficients in light 
of this. 



where v is the frequency in Hertz. As for the spatial filter- 
ing, we apply, in addition to equation (|A.l|h the following: 



G(l) 



1 + exp [I i — I] 1 + exp [I — 1 2] 



(A.2) 



where I is spherical harmonic degree, and Zi , Z2 are the lower 
and upper bounds of the filter, e.g., (h^h) = (0,400) for 
the lowest spatial bandpass of Figure E^a). 
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Appendix A: Filters 

The frequency filter applied to isolate signals in the 5-6 
mHz bandpass is: 



T(y) 



1 



1 



1 + exp 



0.005-M 
0.00005 



1 + exp 



H-0.006 
0.00005 



l,(A.l) 



